THE PLASMONIC EIGENVALUE PROBLEM 



DANIEL GRIESER 

Abstract. A plasmon of a bounded domain Q C R™ is a non-trivial bounded 
harmonic function on R™ \ 8Q which is continuous at 8Q and whose exterior 
and interior normal derivatives at dQ have a constant ratio. We call this ratio 
1 a plasmonic eigenvalue of Q. Plasmons arise in the description of electromag- 

netic waves hitting a metallic particle Q. We investigate these eigenvalues and 
prove that they form a sequence of numbers converging to one. Also, we prove 
regularity of plasmons, derive a variational characterization, and prove a sec- 
ond order perturbation formula. The problem can be reformulated in terms of 
Dirichlet-Neumann operators, and as a side result we derive a formula for the 
shape derivative of these operators. 

Contents 

1. Problem and results 1 
Physics background and previous work 5 
Acknowledgement 6 

2. Basic observations |6 

2.1. Relation to Dirichlet-Neumann operators; proofs of Theorems Q] and [2] 

2.2. Variational principle 

2.3. Scale invariance 

2.4. The two-dimensional case 

3. Perturbation of plasmonic eigenvalues S 
3.1. Analyticity _9 

11 
19 
20 
20 
20 
2C 
21 
21 



3.2. The perturbation calculation 

3.3. Shape derivative of Dirichlet-Neumann operator 

3.4. Splitting at first order 
4. Examples 

4.1. Planar disk 

4.2. Perturbed ball 

4.3. Example: Half space 
References 



1. Problem and results 

Let ft C M™, n > 3 be a bounded domain with smooth boundary dfl and 
connected complement K ra \ O. For functions u_ on Q, and u+ on W 1 \ fi, assumed 
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to be smooth up to the common boundary 9f2, we consider the following problem: 

(1) Au± = in n rcsp. R n \Tt 

(2) u- — u+ = on d£l 

(3) ed n u- + d n u + = on dtt, 

(4) u+(x) =0(\x\ 2 - n ) as |a;|^oo 

where d n u± = n ■ Vw± is the derivative of u± in the direction of the unit normal 
n pointing out of VI. Here e £ M is a constant. Thus, the function u on l n defined 
by U- on VL and by w + on R™ \ f2 is continuous at <9f2, but its normal derivative 
may have a jump, and the ratio of inner and outer normal derivative is required to 
be constant on dtt. It will be convenient to also allow e = oo, in which case ([3]) is 
interpreted as d n u- = 0, d n u + ^ 0. 

We call £ 6 IU {°°} a plasmonic eigenvalue of Q, if the system ((T|)- (J4j) has a 
solution u ^ 0, and then u is called a plasmonic cigenfunction or plasmon for 
£. This name is justified by the origin in the physics of plasmons, and by the 
fact, established below, that e _1 is indeed an eigenvalue of a certain operator on 
d£l. For a given plasmonic eigenvalue we call the linear space of solutions u its 
plasmonic eigenspace. What is usually denoted e in the physics literature is — e 
in our notation. Our choice of sign is motivated by the fact that e > then. A 
corresponding problem may also be formulated for n = 2, but this requires a few 
minor adjustments, see Section [2.41 

Plasmons have been studied extensively in the physics literature, though princi- 
pally for simple geometries, where explicit calculations are possible, and for more 
general domains by numerical methods. See the end of the introduction for more 
on the physics background and references. The purpose of this paper is to put the 
problem in its natural mathematical setting and to bring some more refined math- 
ematical techniques to bear on the study of plasmonic eigenvalues. This will show 
that also for general domains interesting conclusions can be drawn even without 
numerical calculation. In particular, we will employ microlocal analysis to analyze 
the asymptotic behavior of the plasmonic eigenvalues, and differential geometry to 
calculate the perturbation of plasmonic eigenvalues under perturbation of the do- 
main, to second order. This is the order needed for investigating the behavior under 
random perturbations of the domain, a topic of physical interest. Also, the pertur- 
bation results can be useful in obtaining good approximations for non-infinitesimal 
deformations of the sphere, for example, see [3]. 

We now state our results. First, we have basic results on existence and regularity 
of plasmons, and asymptotics of plasmonic eigenvalues. Denote 

V. = {u satisfying ©, ©, ©}, 

and on % consider the scalar product (for simplicity we assume all functions to be 
real- valued) 

(5) (u,v)+ := / Vu • Vw 



where il c = M" \ SI. This is finite since harmonic functions on Jl c satisfying ((4]) 
have Vu(x) = 0(|x| 1_ "). Nonde generacy follows from the unique solvability of the 
Dirichlct problem for f2. Thus H with (•, -) + is a prc-Hilbert space. 
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Theorem 1 (Existence, completeness, asymptotics). Let il C K", n > 3 be a 
bounded domain with smooth boundary. Then the finite plasmonic eigenvalues of 
n form a sequence £1,62,... of positive numbers. For e 7^ 1 the corresponding 
eigenspaces are finite dimensional. 

Furthermore, Ek — > 1 as k — > 00. Also, if E e is the space of plasmons u with 
eigenvalue e then we have an orthogonal decomposition 

e e, 

eSRU{oo} 

Also, £ = 00 is a plasmonic eigenvalue, with one- dimensional eigenspace 

Eoo = {u G H : U- is constant} 



See Section [2T41 for the case n = 2. 

One may also formulate the eigenvalue problem with weaker regularity assump- 
tions. We then have the following result. 

Theorem 2 (Regularity). Let u be a solution of |T]) — ^ with it_ g H s (il), 
u + 6 i/ 1 s oc (ri c ) /or some s > 3/2. 

If e ^ 1 then u is C°° on W 1 \ dQ up to the boundary from each side. 

Here H s is the Sobolev space of order s, and i/j s oc (f2 c ) is the space of functions 
whose restriction to Kr = il c n {\x\ < R} is in H s (Kr) for all i?. The range of s 
is chosen so that the restrictions of u and d n u to the boundary are well-defined by 
the Sobolev restriction theorem. 

Examples show that the statement of the theorem is not true in general for s = 1, 
see Section |4~T1 Note that the only issue is boundary smoothness since harmonic 
functions are always smooth in the interior of the domain. 

The main ingredient in the proof of Theorems Q] and [2] is a reformulation as a 
problem on dfl in terms of Dirichlct-to- Neumann operators, see Section 12.11 This 
is closely related to the fact that for a harmonic function u_ on £1 we have 



(6) / u_d n u_ = / |V U _|^=: Hi 
Jon Ju 

which follows from Green's formula. Similarly, for a harmonic function on M n \fi 
satisfying (j4]), we have 

(7) -/ u+d n u + = I |Vu+| 2 =: + 

Jdn Jn<> 

Theorem [1] was formulated in terms of the scalar product (•, •)+ corresponding to 
the norm || • || + on H. The orthogonality there also holds with respect to the bilinear 
form [u, u)_ := f„ V« • \7v corresponding to j| ■ ||_ which, however, is not positive 
definite - it is zero on E^. 

Equations (J5]) and ([7]) imply that for a plasmon u the plasmonic eigenvalue equals 

(8) e = R p (u), R p (u) : = jjv" + |2 

Denote by H^W 1 ) = {u G Lf oc (W l ) : Vu £ L 2 (R")} the first homogeneous Sobolev 
space. Then we have: 
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Theorem 3 (Variational principle). (|TJ) -((3j) is the Euler- Lagrange equation of the 
functional R p on H 1 {W l ) \ 0. More precisely, the finite plasmonic eigenvalues 
are precisely the positive critical values of R p and the associated plasmons are the 
corresponding critical points of R p . 

Note that zero is also a critical value of R p , with any smooth function with non- 
empty support contained in fi as critical point. Of course, these do not correspond 
to plasmons. 

In this theorem we use H^ oc instead of H 1 since functions satisfying (j4]) need not 
be in L 2 (W l ) for n < 4. 

Next we consider the behavior of plasmonic eigenvalues under perturbation of 
the domain. Thus, fix and let a : d£l -4 1 be a smooth function, and for real 
numbers h close to zero consider the bounded domain tt(h) with smooth boundary 

(9) dQ(h) = {x + ha(x)n(x) : x £ dtt}. 

Thus, O(0) = fl and fl(h) is obtained by shifting the boundary in the normal 
direction by ha. h is the order of magnitude of the shift and a is a 'shape' function. 

In the following theorem, we state the perturbation formulas only in the physi- 
cally interesting case n = 3. Generalization to any n > 2 is straightforward. 

Notation from differential geometry: H and K denote the mean and Gaufi 
curvature of <9f2, Wo = W — HI the trace free part of the Weingarten map W, 
and diva and Va the divergence and gradient operator in the surface dfl. The 
definitions of these quantities arc recalled in Section [3l Also, 

(10) (g,g') := / 99' dS, g,g'eC°°(dn) 

where dS is surface measure (all functions are assumed real- valued) , and for an 
operator A on C°°(d£l) we let 

(g\A\g'):= (g,Ag') 

Similarly, for vector fields V, V on dfl and operator A acting on vector fields, 
(V\A\V) := (V, AV) := J gn V ■ AV'dS. All operators appearing below are 
self-adjoint. 

Theorem 4. Let e ^ 1 be a finite plasmonic eigenvalue of £1 C K 3 with eigenspace 
E. Then there are ho > and real analytic functions h i— > eW(ft), h i— > u^\h) 
defined for \h\ < h , i = 1, . . . ,dim_E, such that for each h the numbers e^(h), 
i = 1, . . . , dimi? are plasmonic eigenvalues ofil(h) with eigenf unctions u^ l \h), and 
uW(0),...,u( dim£ )(0) are a basis ofE. 

For a fixed analytic branch s{h), u(h) of eigenvalue and eigenfunction with e(0) = 
e and \\u(h)\\- = 1 for each h we have for i :— ^e(O), ii(x) := j-^u(x,0), e := 
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\dh) 
(11) 



(12) 



:(0) the formulas: 



qi(u) := (e + 1)[ - (Vau_ | a | V 9 w_ ) + e(d n U- | a | <9„it_ )] 
(e+1) / a[-\V a u_\ 2 + e{d n u_) 2 } 



i on 

>- 

[d n u- 



-ia - (e + l)a 2 W 
(e 2 - l)a 
-(e + l)a 

— £ 

ae(e + (e + l)aff) 
e(e + l)a 



V a (a<9„u_)) 

a„?i_ } 

d n u- ) 



Here, u is a solution of the inhomogeneous system 
Aii = 

(13) ?i- — ii + = — (e + l)ai9„u_ 

ed n ii- + d n u+ = —ed n u- + (e + 1) diva(aVaw_) 

^4feo, m(0) is a critical point of the functional R' p {u) 
in wit/i critical value i. 



II 



m M 3 \ Sf2, 

on dQ 
on E, with q\ defined 



The meaning of analyticity of h \ 



(h) will be explained in Section |3~T1 Note 



that m in (|13p is not unique, but the value of e calculated from tjl 2|) is independent 
of the choice of iz. There is a similar inhomogeneous system determining ii, see 
(f2Tj) -([23 |) . In the proof of the theorem, the formulas for e, e are obtained from the 
requirement that these systems have a solution. The last statement of the theorem 
helps to identify which u to use in (TlTj) in case of degeneracy, i.e. dimi? > 1. 

As a side result we also obtain a perturbation formula for the Dirichlet-Neumann 
operator which may be of independent interest, see Theorem ll3l This is sometimes 
called the shape derivative of the Dirichlet-Neumann operator. Recently there has 
been some interest in such results in the context of water waves, see [llj for example. 

Physics background and previous work. The problem QJ-lHl) with n = 3 
arises in the physics of surface plasmon polaritons. These are waves which couple 
electromagnetic fields to the electron gas of a metal, and thus propagate along the 
interface of a metal and a dielectric material or vacuum. As outlined in several 
recent reviews, the particular interest in these excitations stems from the fact that 
by altering the structure of the metal's surface, the properties of these plasmons, 
including their interaction with light, can be specifically tailored [Tl [2l [T2l [TBI 114] . 
This offers a high potential for developing new types of photonic devices, combining 
the speed of photonics with the nanoscale dimensions of electronics. Possible ap- 
plications include, among others, plasmonic chips, data storage, nanolithography, 
and highly sensitive molecular detectors. 

Our setting models a metallic particle whose shape is given by f2, surrounded 
by vacuum. The function u is the potential of the electric field, and the equations 
arise from Maxwell's equations in the quasistatic approximation, which assumes 
constancy of the magnetic field, and is justified if the wavelength is much larger 
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than the dimensions of the particle, which is the case for particles on the nano scale. 
See [El Chapter 7] for details. In the context of surface plasmons the quasistatic 
approximation has often been used in the literature, see [HI |5] for example, and its 
validity has been verified experimentally. 

The boundary conditions ([2]), (|3]) arc obtained from continuity of the tangen- 
tial component of the electric field and of the normal component of the electric 
displacement, see Chapter 4.4]. The number — e is the relative permittivity of 
the material of the particle, which is coupled to the frequency by a law specific 
to the material. Thus, plasmonic eigenvalues correspond to frequencies of elec- 
tromagnetic waves at which there is a strong excitation of the electron plasma in 
the metallic particle fl. The name plasmonic eigenvalue was coined in [5], where 
also the asymptotics Ek — > 1 was proved. The paper [3] contains the first order 
perturbation formula (jlip . 

The spectrum of the Laplacian on a bounded domain, and its relation to the 
geometry of the domain, has received been studied intensively for many years, so 
it may be worth while to compare it to our problem. There are similarities, for 
example the existence of discrete spectrum and the presence of the low energy 
and the high energy regime, which require very different mathematical techniques. 
There are also clear differences, for example the asymptotic behavior and the fact 
that plasmonic eigenvalues are scale invariant while the Laplacian eigenvalues are 
not. 

Acknowledgement. I am grateful to Gohar Harutyunyan for checking the per- 
turbation formula (|12[) independently. 

2. Basic observations 

2.1. Relation to Dirichlet-Neumann operators; proofs of Theorems [l] and 

[2j Theorems [T] and [2] are proved by reformulating the plasmonic eigenvalue problem 
as an 'honest' eigenvalue problem for an operator on C°°(dSi) which is built from 
Dirichlct-to-Neumann operators. For basics on harmonic functions in Q and R n \f2, 
and for proofs of some of the facts stated below, see for example [4] . The behavior 
at infinity of harmonic functions may, for example, be analyzed by using a Fourier 
decomposition on large circles. 

The (interior) Dirichlet-to-Ncumann operator of f2 maps Dirichlct to Neumann 
data of a harmonic function. That is, 

TV- : C°°(dn) -> C°°(dn), g i-> «9„u_ 

where w_ is the unique solution of the Dirichlct problem Aw_ = in £1, U-\qq = g. 
Similarly, the exterior Dirichlet-to-Ncumann operator is 

Af+ : C°°{dn) -> C°°{dn), g^d n u+ 

where it+ is the unique solution of the Dirichlct problem Au+ = in R n \il, u + \qq = 
g, u+(x) = 0(\x\ 2 ~ n ) as \x\ — > oo. It is standard that this decay requirement 
already follows from boundedness of u+. 

JV± are self-adjoint operators on L 2 {dU) with the standard scalar product (fTO)) . 
The symmetry of Af- follows immediately from Green's formula, and for M+ can 
be reduced to that of Af- by use of the Kelvin transform. 

Note that while Af- is a positive operator by (JB]), Af+ is negative because the 
normal n points into K™ \ fl rather than out. 



THE PLASMONIC EIGENVALUE PROBLEM 



7 



Proposition 5. The plasmonic eigenvalue problem (JTJ) is equivalent to the 
generalized eigenvalue problem 

(14) (eA/1 + N+)g = 0, 

in the sense that smooth solutions (e,u) of (jTJ) -((4j) are in 1-1 correspondence with 
smooth solutions (e,g) of (|14j) . via g = u\qq. 

Proof. By the previous discussion we have a 1-1 correspondance 

n <— > c°°(dn) 

Under this correspondence ([3]) is equivalent to (fT4)) . □ 
For e C°°(<9fi) we denote 

(g,g')+ :=-{gM) 

If u,u' € Ji are the functions on R n corresponding to g, g' under the correspondence 
(US]) then Green's formula, cf. (0, gives 

(15) (g,g')+ = (u,u') + 

where the latter is defined in ([5]). This justifies the notation and shows that (-,•)+ 
is a scalar product on C°°(d£l). The latter fact also follows from the well-known 
fact that Af+ is invertible. 

The problem (|14p can easily be transformed into a standard eigenvalue problem: 
The plasmonic eigenvalues are the reciprocals of the eigenvalues of the operator 

A = -N+ 1 N- 

This operator is symmetric with respect to (•,•)+• Introduce the Hubert space 
H := completion of C°°(90) with respect to (•,•)+■ 

Proof of Theorem^ Positivity of plasmonic eigenvalues follows from ([5]): R p (u) 
must be non-zero since otherwise w+ = 0, so u_\qq = 0, hence u = 0. 

We use some concepts and facts from microlocal analysis. See [4] for background. 
It is well-known that Af-,JV+ are first order classical pseudodifferential operators 
on C°°(<9f2) with principal symbols 

<7(Ar-) = Kl> ^ + ) = -iei. 

The fact that A/+ is an invertible elliptic first order operator implies that % = 
H^idfl), the Sobolev space of order 1/2. By the standard theory of pseudodiffer- 
ential operators, A = — Af+ 1 Af- is a zeroth order pseudodifferential operator with 
principal symbol — (— 1^| _1 ) • |£| = 1, hence A = I + R where R is a pseudodifferen- 
tial operator of order — 1. Therefore, R is a compact operator on TL, and since A is 
sclfadjoint, so is R, so R has real nonzero eigenvalues ri, T2, • • • — > of finite multi- 
plicity, and may have a zero eigenvalue r$ = of any multiplicity. The eigenspaces 
span H and are orthogonal with respect to (•,■)+■ Then ak = 1 + ru are the eigen- 
values of A, so a^T 1 are the plasmonic eigenvalues and a^ 1 — > 1 as k — > oo. By (p~5]) 
the plasmonic eigenspaces are orthogonal with respect to the scalar product (•, •)+ 
on H.. There is an eigenvalue = which corresponds to the kernel of 7V_, that 
is, the constants in Ji, and these correspond to the space E^. □ 



8 



DANIEL GRIESER 



Proof of Theorem Clearly, the correspondence in Proposition [5] also holds for 
.ff s -plasmons. The operator eAf- +7V+ has principal symbol (e — 1)|£|. If e ^ 1 
then this is invertiblc for each £ ^ 0, hence the operator eJ\f- +Af+ is invcrtible then. 
So by elliptic regularity its kernel consists of smooth functions. These functions are 
the boundary values of the corresponding plasmons, so these are smooth up to the 
boundary by elliptic boundary regularity. □ 

2.2. Variational principle. 

Proof of Theorem^ We calculate the first variation of R p . Let u 6 iJ 1 (R") be 
such that Jq |Vm| 2 > and let e = R p (u). For v <E Cp°(R n ), the space of smooth, 
compactly supported functions, we have 

(16) ^ (v) = 2 ( [ Vu + Vv + - e [ Vu. 

(17) 



/ n |V«-| 



( / (d n u + + ed n u-)v - I (Au+)v+ + e / (Au-)v- 



using integration by parts. If u satisfies (HJ-© then it follows that = for 

all v G Cg°(R n ). Now C£°(IR' 1 ) is dense in H x (It n ), and for it satisfying in addition 
(HJ) the right hand side of (JTBJl depends continuously on u, so m is a critical point 
of R p , and R p (u) = e > by definition and Theorem [T] 

Conversely, assume w is a critical point of R p with e = R p (u) > 0. Using smooth 
test functions v localized near a point in f2 or the interior of fi c one sees from (|16p. 
(fT7|l that Au = on K™ \ and then using v localized near points on <9f2 one gets 
ed n U- + d n u + = at dfl. By the Sobolcv restriction theorem we have u_ = u + at 
dVl. Finally, a harmonic function u on Q, c with Vu G L 2 (f2 c ) satisfies ((3]), so u is a 
plasmon with eigenvalue e. □ 

2.3. Scale invariance. As opposed to the eigenvalue problem for the Laplacian, 
the plasmonic eigenvalue problem is scale invariant: If t > then the domains fl 
and tQ = {tx : x G f2} have the same plasmonic eigenvalues, since an cigenfunction 
u for CI turns into an cigenfunction ut(x) := u{x/t) for tQ with the same e. 

2.4. The two-dimensional case. The condition (j4} on the behavior of u at in- 
finity is natural for all dimensions n > 2 in that it gives unique solvability of the 
Dirichlet problem on fl c for any boundary data on dfl. However, this condition 
is not suitable for formulating an analogue of the plasmonic eigenvalue problem 
for ?? = 2 because any constant function u will satisfy (JTJ)- f° r an U value of e. 
On the other hand, adding a constant to a plasmon u with arbitrary plasmonic 
eigenvalue will yield a plasmon with the same eigenvalue. Therefore, it is natural 
to consider plasmons as elements of the quotient space H/ {constants}. Since any 
bounded harmonic function on il c C M 2 is of the form constant + 0(|x| _1 ), this is 
equivalent to replace the decay condition ((4|) by the condition 

(0k=2) u + {x) = 0{\x\^) as \x\ ->• oo 

Denote the space of plasmons satisfying ([1]), and (|4I-9j) by %' . On this space 
(-,-)+ is a scalar product, and Theorem Q] holds, except that now we have the 
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decomposition 



and there is no space E^,. Theorems [2] and [3] hold as in n > 3, with the same 
proofs. 

In Section 12.11 the definition and discussion of the Dirichlet-Neumann operators 
jV± before Proposition [5j and also Proposition [5j hold for n = 2 (where condition 
(H|) is used to define Af+). However, the operator Af+ is not invertiblc if n = 2. 
More precisely, we have 



and by self-adjointness both ranges are equal to the orthogonal complement of the 
constants (with respect to the standard scalar product ([TP])), which we denote by TZ. 
Therefore, both A/+ and A/1 descend to invertible operators on C°°(dfl) /{constants}, 
which may be identified with TZ. Then the plasmonic eigenvalues are again the re- 



C°°(<9f2)/ {constants}. The operator A is invertible, hence oo is not a plasmonic 
eigenvalue. The proof of Thcorcm[T]carrics over if one extends A to % by continuity 
and by setting it zero on constants. 

There is a subtlety, which is the reason why quotient spaces are more natural in 
this discussion than their replacements %' and 1Z: The boundary values of harmonic 
functions on fl c satisfying condition (|4ln — arc not the functions in 1Z. Rather, they 
are characterized as follows. There is a unique harmonic function v on Q c having 
boundary values v\qq = and satisfying v(x) ~ ^- log \x\ as \x\ — > oo. This function 
may be constructed by choosing any point y € fl and letting v = G — G' where 
G(x) = y~ log \x — y\ is the Newton potential centered at y and G' is the unique 
bounded harmonic function on tt c having the same boundary values on <9f2 as G. 
Let g = d v v e C°°(dD,). Then we have: 

g G C°°(90) is the boundary value of a harmonic function on f2 c 

satisfying (|4l,-9|) if and only if (g,go) = 0. 

Necessity follows from Green's formula, and then sufficiency from a dimension ar- 
gument. The function go is not constant, except when f2 is a disk. This means that 
the harmonic extension u of g € 7Z will not satisfy the bound (|4L — ?p in general. In 
order to obtain a plasmon satisfying (|4l, -?|) one has to add a constant to u. 



In this section wc prove Theorem 01 

3.1. Analyticity. The analytic dependence of the plasmonic eigenvalues on h is 
proved by standard methods: First, translate into a problem on a fixed domain, 
then show that the operators involved depend analytically on h, then apply standard 
perturbation theory. 

Wc provide the details: In a first step, we translate the problem for fi(/i) to a 
(/i-dcpendent) problem for the ft,-indcpcndcnt domain tt. To do this, we extend the 
vector field an on <9S1 to a smooth vector field V on R™, which vanishes outside 
some small neighborhood U of dSl. This can be done by first extending an to a 
neighborhood of dfl and then multiplying by a smooth function on R" supported in 



kerA/+ = ker Af- = {constants} 




3. Perturbation of plasmonic eigenvalues 
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this neighborhood and equal to one on dtt. Then define maps Th : K™ —> M. n ,x i-> 
x + /iV(a;). For any h, Th is the identity outside U, and by compactness of d£l and 
the inverse function theorem it is a diffcomorphism on R" for h sufficiently close to 
zero. By definition, Th is a diffcomorphism 51 — > f2(/i). 

Now the map 7), is used to 'transport' data related to Q(h) back to J7: Let g cuc i 
be the Euclidean metric on R" and gh = T h g cuc i its pull-back. Denote by Ah, V/, 
the Laplacc-Bcltrami and gradient operators for the Riemannian metric gh on R" 
and by the normal vector field with respect to gh on <9f2. Then the plasmonic 
eigenvalue problem for Cl(h) is equivalent to the problem for Q with A, d n replaced 
by Ah, d nh . More precisely, u solves the former if and only if T£u := uo Th solves 
the latter. Also, in the discussion of the variational principle one needs to replace V 
by V/i, and in Proposition [5] the Dirichlct-Neumann operators need to be replaced 
by those for the metric gh, denoted A/±,fe. 

We can now state the precise meaning of analyticity in Theorem 2) 

Proposition 6. Let e =/= 1 be a finite plasmonic eigenvalue of Q with eigenspace 
E . Then there are ho > and real analytic functions 

e^' : (—ho, ho) — > R, u^> : (—ho, ho) -> {u e C(R n ) : and u\ Rn \Q are smooth} 

for i = 1, . . . , dimi? such that ^'(O), . . . , u^ dlmE '(0) form a basis of E and the pairs 
((T^VuW^.eW^)) satisfy the pi asmonic eigenvalue problem for the domain 
Q(h), for each h <E (—ho, ho). 

The analyticity of uS 1 ' in the given space is to be understood as analyticity of 
the restrictions in the Hilbcrt spaces H s (£l) and H s (£l c C\ {\x\ < R}) for each R > 
and s e R. 

Proof. We first consider the reduction of the problem given in and after Proposi- 
tion [5] Thus, we are dealing with an eigenvalue problem for the operator Ah = 
—J\f7 1 h Af-.h on L 2 (dfl). The analytic dependence of plasmonic eigenvalues and 
eigenfunctions, as elements of L 2 (dfl), follows from standard perturbation theory, 
see [7], if we can show that h i-> A/± h are analytic families of operators. Wc 
show this for Af h := A/L^, the argument for is analogous. By [7J Theo- 

rem VII. 4. 2] we only need to show that the associated family of quadratic forms, 
= (9,Nh9), is analytic in h, in the sense that the form domain is inde- 
pendent of h and that h — > qh(g) is analytic for each g in this domain. Now 
since Afh is an elliptic pscudodiffcrcntial operator of order one, the domain of q h 
is the Sobolev space H l / 2 (d£l) for each h. Also, for each g e H x / 2 (dVt) we have 
(compare (O) (g,Nh9) = Jo \^h^h\ 2 dvol^ where Vh € iJ 1 (f2) is the solution of 
the Dirichlet problem on Q with boundary value g (here we need to impose also 
the decay condition in the case of Af + .h). So we need to prove analyticity of 
h i v (g,J\Th9) = Jo l^/i^Til 2 dvol^. We first prove analyticity of h v h : Since 
the metric gh depends analytically on h (in fact, Sh( x ) is given by the matrix 
(I + hDV(x)Y(I + hDV(x)) in standard coordinates), so does A^. Let Gh be 
the inverse of A^ on f2 with Dirichlet boundary conditions at dfl. Then Gh is 
analytic in h by [3 Theorem VII. 1.3]. We have = v — GhAhV where v = i>o is 
the unperturbed solution. This shows that Vh is analytic in h. Finally, V ' hVh and 
dvol/j = \J dot gh dx are analytic in h, so we have proved the analytic dependence 
of Nh on h. 
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Now the cited theorems give the existence of (h) and of branches g^ (h) <G 
L 2 (dfl), depending analytically on h. Then since eW(0) = e ^ 1 we can choose 
ho small enough so that £«(ft) ^ 1 for all h £ ( —ho, ho). Then the standard 
parametrix construction gives, for any s G N, a left parametrix of P/j = 
e {i) {h)N-, h + M+.h with error term P h = Q ft P ft - 7 G ^- s {dQ), with Q h de- 
pending on h analytically. Then = QhPhg^(h) = g^'(h) + RhgW(h), and this 
shows that g^ l \h) is analytic with values in H s (d£l). Finally, this implies that 
u^(h), defined as the solution of the Dirichlet problem with respect to with 
boundary values g( l '(h), depends analytically on h in any Sobolev space. □ 

3.2. The perturbation calculation. 

3.2.1. Procedure. Wc first explain our procedure. Consider ([I])-(j4]) with all quanti- 
ties depending on h. Differentiate these equations in h and set h = 0. This yields 
a system of equations 

(18) Ail = in M 3 \ dtt 

(19) ii^-u + = F on 9fi 

(20) ed n u- + d n ii + = —id n u- + G on d£l 

satisfied by ii, where F,G are explicit expressions in u and e, see (|4Tj) . (|42]) . This 
is an inhomogeneous version of the system (HJ)— - Since the homogeneous system 
has a nontrivial solution (the given u), solvability of the inhomogeneous system 
implies certain compatibility conditions for the right hand sides F, —id n + G, see 
Lemma [T2l This yields a formula for e. Similarly, differentiating twice in h before 
setting h = gives 

(21) Ail = in M 3 \ dfl 

(22) it_ - u + = F 2 on dft 

(23) ed n u- + d n u + = —ed n U- + G% on dft 

with i*2, G2 given in (|43|). (|44|) . and in the same way as before this yields the formula 
for e. 

3.2.2. Calculation of derivatives at <90. In order to determine F, G, F2, Gi in (|19|) - 
(|23|) . we need to calculate certain derivatives at the boundary. It is important to 
express these in invariant differential geometric terms. For a point p € dfl we 
decompose orthogonally M 3 = T p M. 3 = T p dfl span n(p) and express all quantities 
with respect to this decomposition. For example, for a function u we have 

Vu = Vqu + (d n u)n. 

Denote by n(q, h) the unit normal to dSlh at q G d£lh- 

Proposition 7. Fix p € d£l. Denote by n{p), n(p) the first and second derivative 
in h of the unit normal n(qh, h) at = p + ha(p)n(p) E d£lh- Then we have at p 

h = — Vga 

h = -W{V d a 2 ) - \V d a\ 2 n 
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Furthermore, if x i— > u(x, h) is harmonic on flh for each h, smooth in (x, h) up to 
the boundary, then h t— > Vu(qh, h) has first and second derivatives at h = 

(Vu)' = a(T\u + T[d n u) + a(N x u + N[d n u) n +V« 

(Vtt)" = a 2 (T 2 u + T^d n u) + a 2 (N 2 u + N^d n u) n + 2(Vu)'-Vu 

where 

Tx = WV d> T[ = Va, 

Ni = -A 9 , N[ = 2H, 

N 2 = 2(-HA d - div a WV 9 + V d H • V a ), ^ = - A fl + 8i? 2 - 2A'. 

Here X is the Gauss curvature of dfl and Ag = div^ Va is the Laplace-Beltrami 
operator on dQ. The definitions of these and the other quantities is recalled in the 
proof of Lemma [5] The operators TJ, T/, iVj, N[ act on functions defined on d£l. 

The tangential component of (Vu) is not needed, so we don't calculate it. The 
term (Vu)' can be calculated using the same formula as for (Vu)' since ii(h) is also 
harmonic for each h. The same conclusion holds if «(•, h) is harmonic on R™ \ flh 
for each h. 

The main point in the formulas for the variation of Vu is that they involve 
at most first derivatives in the normal direction. The second and third normal 
derivatives, which occur in calculating the left hand side, are expressed in terms of 
tangential derivatives and curvature terms, using the harmonicity of u. See Lemma 

ma 

Proof. Fix p <G d£l. We introduce coordinates in K 3 so that p = and that the 
tangent plane to dfl in p, T p dfl, is horizontal and f2 lies below it. Then we can 
parametrize the surface dfl near p as a graph: 

(24) x(v, w) = (v, w, f(v, w)), 

so that, near p, dSl = {x(v,w) : v,w near 0}, and / = 0, V/ = at (0,0). We 
may rotate coordinates so that / has the form 

(25) / ( WiU ,) = ^ t ,2 + | U7 2 +0 3 

for some A,/j£i, where here und in the sequel O k denotes any quantity vanishing 
at least to order k at v = w = 0. Geometrically, A,/i are the principal curvatures 
of dfl at p. 

Let n = n(v, w) be the upward unit normal: 

(26) n(v, w) = (-/>.™W»fr."0.1) 

y/1 + f*(v,w) + f*(v,w) 

Expanding this to second order we get 

n(0,0) = (0,0,1), 

(27) n„(0,0) = (-A,0,0), 

^ui;(0, 0) ( fvvvi fvvwi A ). 

and analogously with v, w interchanged and A replaced by /z. Here and in the sequel 
derivatives of / are evaluated at (0,0) if no arguments are given. 
We now express geometric quantities of dfl in the coordinates v, w. 
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Lemma 8. For the metric tensor g, second fundamental form II and Weingarten 
map W we have in the coordinates v, w 

(28) g = Id + O 2 

(29) n = w + o 2 = ( * ) + ( v { mw + wf ; vw v { vvw + w { ututu ) + O 2 

y U /i / y VJ VVW -\- WJ VWW VJ WWV t WJ WWW J 

Here Id is the identity matrix. The derivatives of f are evaluated at the origin. 
The intrinsic differential operators in the surface are 

(30) V 9 = V + O 2 , diva = div +0\ A d = A + O 1 
where V,div, A are the standard (Euclidean) operators in v,w. 

We note the following consequences of the lemma, needed below: The mean 
curvature is 

H = i traced = -(A + /i + vAf v + wAf w ) + O 2 . 
and then simple calculations give 

(31) diva WV e = Xd 2 + M d 2 + V(A/) ■ V + O 1 
and 

(32) V d H = ^VAf + O 1 . 

Proof. For the metric tensor g we have 

gn gi2 \ _ ( d v x ■ d v x d v x ■ d w x 
92i 922 J V d w x ■ d v x d w x ■ d w x 

1 + fl{v,w) f v (v,w)f w (v,w) 
f v (v,w)f w (v,w) l + f 2 (v,w) 

The second fundamental form, expressing curvature, is defined by 

j j _ ( H\\ H\2 \ _f n(v,w) ■ d 2 x(v, w) n(v,w) ■ d 2 w x(v,w) 
\ Ihi H22 ) \ n(v,w) ■ d 2 w x(v,w) n(v, w) ■ dl,x(v, w) 

where 

/ \ o2 / \ (~fv(v, w),-f w (v,w),l) . . fvv(v,w) 

n{v,w) ■d;x(v,w) = , = -. = =■ (0,0,f vv (v,w)) - 



g{v,w) 



= Id + O 2 



^l + f 2 (v,w)+f 2 (v,w) y/1 + f 2 + Jl 

= f vv (v,w) + 2 , 
n(v,w) ■ d 2 w x(v,w) = f vw (v,w)+<D 2 , 
n(v, w) ■ d^x(v, w) = f ww (v,w) + O 2 . 

Expanding to first order in Taylor series gives the formula for // in (|2"9")l . 

The Weingarten map (or shape operator), again expressing curvature, has com- 
ponents 

2 

W i =Y,Hik9 k3 i,j I--'. 
fc=i 



where ( g k i ) is the inverse matrix of p. From (|28]l it follows that g 1 = Id+O 

V / k.j — 1,2 

and hence W = II ■ g~ l = II + O 2 , so d29j is proved. 
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The gradient of u in the surface is 

v,,=(»:: £ ){z)={^){:: ) + & - 

Next, g = Id + O 2 implies ^/det g = 1 + 2 , so for a vector field Z = (Z V ,Z W ) 
the divergence in the surface is 

div 9 Z = (d v (y/dct gZ v ) + d w ( ^dct gZ w )^j 

= d v Z v + d w Z w + O , 
and this gives also the result for the Laplace-Beltrami operator, Ag = divg Vg. □ 

We now prove the formulas for n,n in Proposition [Jj We parametrize d£lh by 

y(u, iu, ft,) = x(y, w) + ha(v, w)n(v, w) 

with x(v, w) as in (|23]l. For p = 0we have q h = F(0, 0, ft), and iV(ft) = Y" t ,(0, 0, ft) x 
Y w (0, 0, ft) is a non-normalized normal to dVlh at this point. From (f27j) we have 
1^,(0,0, ft) = (l — a\h,0,ha v ) and analogously 1^,(0,0, ft) = (0, 1 — afj,h,ha w ), so 
(with a and its derivatives taken at (0,0)) 

N(h) = (—a v h + aa v /ih 2 , —a w h + aa w Xh 2 , 1 — a(/x + A)ft + a Xfih ) 
= A + Bh + Cti 2 

where A — (0,0,1), B = (—a v ,—a w ,—2Ha), C = (aa v [i,aa w X,a 2 K). Expanding 
n(qh, ft) = pm^Ti in Taylor series results in n(qi l , ft) = 

A+(b-A(A-B)) h+ (c - B(A ■ B) - A^f- +A-C-^(A- 5) 2 ) j ft 2 +0(ft 3 ). 

The coefficient of ft is (— a v , —a w ,0) = — Vga by (|30[) . and the coefficient of ft 2 /2 is 
(-2aa v X, -2aa w Lt, -|Va| 2 ) = -W(\7 d a 2 )- |V 9 a| 2 n by d29j) and (|30l) , which proves 
the first part of Proposition [7J 

To prove the second part of Proposition [JJ it is convenient to use coordinates 
normal with respect to dfl, near p = 0. That is, we parametrize a point X = 
(Xi, X2, X3) near p by coordinates (w, w, s) where s is the signed distance from X 
to dil and x(w,u>) is the point on <9f2 closest to X. That is, 

X(v, w, s) = x(i>, w) + sn(v, w). 

Lemma 9. The derivatives in the coordinates X\, X2, X3 and v,w,s at the point 
p = are related as follows: 

dxx = d v , dx 2 = d w , dx 3 = d s ; 
d 2 Xl =d 2 v -Xd Sl d% 2 =dl-fxd s , d 2 X3 =d 2 ; 
(33) d Xl dx 3 = d v d s + Xd v , dx 2 dx 3 = d w d s + ^id w ; 

d Xl dx 3 = d 2 d s - Xd 2 + 2Xd 2 - X 2 d s + f vvv ■ d v + f vvw ■ d w ; 
d x 2 dx 3 = d 2 ,d s - \id 2 s + 2\id 2 w - y 2 d s + f www ■ d w + f wwv ■ d v . 
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Proof. By differentiating and using the chain ruie repeatediy, one gets 
K ") x =n-d x , d 2 s = n ■ n ■ d 2 x ; 



dx 

d v = X v ■ dx, d 2 — d v (X v ■ dx) = X v ■ X v ■ d x + X vv ■ dx ; 
d v d s = n ■ X v ■ d 2 x + n v ■ d x , 

d 2 d s = n ■ X v ■ X v ■ d x + n ■ X vv ■ d\ + 2n„ • X v ■ d\ + n vv ■ d x , 

and analogous formulas for the w-derivative^]. 

We now evaluate all quantities at p = 0, i.e. v = w = s = 0. Using (p?T|) and 
X v (0, 0, 0) = (1, 0, 0), X vv (0, 0, 0) = (0, 0, A) we get at p 

d a = dx 3 , d 2 s = d 2 x 3 ; 

d v = d Xl , d 2 = d Xl + \d Xs ; 

d v ds = d Xl dx 3 - \d Xl , 

d 2 d s = d Xl dx 3 + ^9x 3 — 2Ad 2 fl — fvvvdx t — fw W dx 2 ~ A 2 <9x 3 - 

Analogous formulas hold with v replaced by w, A by fx and with dxi , dx 2 in- 
terchanged. Solving these equations for the derivatives dx t one gets the desired 
relations. □ 

Lemma 10. Let u be harmonic on Q, smooth up to d£l. Then 

(34) (n-V) 2 u = N lU + N[d nUl 

(35) (n-Vfu = N 2 u + N' 2 d n u 

at p = 0, with Ni, N[ as given in Proposition^ 

Here \(n ■ V) 2 it] (p) is understood as (n(p)V) 2 u, evaluated at p; that is, n is 
not differentiated^- In the coordinates above, we have at p = 

(n • V) 2 u = d X3 u, {n-Vfu = d X3 u 
Proof. A short calculation, using (|53")h ([30]) and d s = dx 3 = d n , gives, at 0, 
dl 3 u = -{d 2 Xl +d 2 X2 )u 

= -{dl-\d s )u-(d 2 w - ^d s )u 
= -(d 2 + dl)u+(X + fi)d s u 
= —Agu + 2Hd n u. 
Next we evaluate d Xa u at 0. By the last formulas in ([3"3"| we get 
d 3 x 3 u = -{d 2 Xl +d 2 X2 )dx 3 u 

= -(d 2 d s + d 2 w d s )u + (\ + f i)d 2 u 

- 2{\dl + iid 2 w )u+{\ 2 + pL 2 )d s u 

\fwv®v H~ fwww&w H~ fvvw&w H~ fwivv^v)^- 
V(A/)-V 



3 

For vectors V, W we write V ■ W ■ d x = Yl ViWjdx dx , for short, and similar for third 

M=l * 3 

derivatives in X. 

2 Of course it would not make sense to differentiate it since it is only defined at 8Q. But this 
is in contrast to the consideration of (d n u)' below, where the normal depends on h, hence must 
be differentiated. 
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Using (00]), d31]), j32j and then <9 S = d n (everywhere), d 2 s u = d X3 u = (-A 9 + 
2Hd n )u (at p only) we obtain 

d 3 Xs u = -A a d s u + 2Hd 2 a u + (4H 2 - 2K)d s u - 2div a (W(V a u)) + 2V a H ■ V a u 

= ( - A a + 8H 2 - 2K^d n u + 2^- HA a - div a (W(V a )) + V a # • Vd)u 

as claimed. □ 

Now we prove the second part of Proposition [7J We have 4r [(Vu)(qh,h)} = 
(qii ■ V)(Vu)(q h ,h) + (Vu)(qh, h). With q h =p + ha{p)n(p) we get q h = a(p)n(p). 
At p = 0, n(p) ■ V = <9av So we get, using (|33]l and (|34|) . with everything evaluated 
at p = 0, h = 0, 

dx 3 d Xl u \ 

(Vw) = a(n • V) Vu + V« = a I dx 1 dx 2 u +V« 

^« / 

/ \ 1 

Vu 

The hrst term is a\7 a d n u by <j30j), and the second is aW(\7 a u) by <(29j) , (|30|) . 
Finally, 

(Vu)" = a 2 (n- V) 2 Vw + 2a(n • V)Vu + Vu 

and the normal component of the first term is a (d x u)n = a 2 (N 2 u + N^d n u)n by 
(|35|) while the last two terms can be rewritten as 2(Vit)' — Vu. This hnishes the 
proof of Proposition [71 □ 

Lemma 11. Let u be as in Proposition^ The derivatives ofh^ u(qh,h) at h = 
are 

(36) (it)' = ad n u + ii, (u)" = —a 2 A a u + 2a 2 Hd„u + 2ad n ii + it 
The derivatives of 7i i— > (d n u)(qh, h) at h = are 

(37) (S n u) ' = Pi («) + Qi (SU) + 9„ ii 

(38) (cU)" = P 2 (u) + Q 2 (d n u) + 2Pi(u) + 2Q 1 (d n ii) + d n u 
with 

Pi = -div s aV 9 , 
Qi = 2aH, 

P 2 = -2(a 2 HA a + div d (a 2 W)V d -a 2 (V a H)-V t 
Q 2 = -div s a 2 V s + a 2 (8P 2 -2K)-\V a a\ 2 : 



Proof. The first equation of (|36|) is obvious. The second one follows from (u(qh, h)) = 
a 2 (n ■ V) 2 u + 2a(n • V)u + u and (j34|) . 

Equation ((37)) follows from (d n u) = (n(g^,/i) -Vu(qh, h)) = n-Vu + n- (Vu) , 
and Proposition [7j The terms acting on it are 

-(Vgo) • V a + aA^i = -(Vao) ■ V a - aA s = - div a aV a = Pi 

and there is one term acting on d n u, which is aN[ = Q\. 
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For (|38]l we write (dr.it) = (n(qh,h)-S7u(qh,h)) = ri-Vw+2n-(Vw) +n-(Vu) , 
use Proposition [7] and sort by terms involving u, d n u, u, d n u. The terms acting on 
u are 

- W{V g a 2 ) ■ V a - 2{V d a)aT 1 + a 2 N 2 
= - (V a a 2 ) • WVg - (V 9 a 2 ) • WV g + a 2 N 2 (since W is selfadjoint) 
= - 2 [(V 9 a 2 ) • M^V 9 + a 2 div 9 iyV 9 + a 2 HAg - a 2 (V 9 ff)V a ] = P2 
The terms acting on d n u are 

- |V s a| 2 - 2(V a a)aT[ + a 2 N' 2 
= - |V 9 a| 2 - (V a a 2 ) • Vo - a 2 A 9 + a 2 (8H 2 - 2K) = Q 2 

The terms involving it, and d n ii arise from one of the ft-derivatives in (n(qh,h) ■ 
X/u{qh, hy\ falling on the second argument of u, hence arc twice the corresponding 
terms in (d n u)' . □ 

3.2.3. End of proof of the perturbation formulas. Now let u(-,h), e(h) be a solu- 
tion of (|TJ)-(|4j) for f2/j, for each ft, smoothly depending on all variables (up to the 
boundary of £lh from either side). 

We first prove iP] ) -([25] ) and calculate F, G, F 2 , G 2 . If x <£ dVL then x <£ dVL h for 
sufficiently small ft, and (fl8|) . (|2"Tj) follow. The boundary conditions are 



(39) u-(q h ,h)-u + (q h ,h) = 

(40) e{h)d n u_(q h , ft) + d n u + {q h , ft) = 

where qh = P + ha(p)n(p), for each p € <9f2 and ft. We differentiate (f39]) and 
use the first equation in (|36p for u_ and u+, at ft, = 0. This gives u_ — u + = 
—a(d n u- — d n u + ). Using (|40[) at ft = we obtain (fT9|) with 

(41) F = -(e+ l)ad n u_ 

Next, we differentiate (|40p at /i = and use ||3"T|) for u_ and u+. We obtain 

= id n u- + e(d n u-)' + {d n u + \ 
= id n u- + P\(eu- + u + ) + Qi{ed n u- + d n u + ) + ed n U- + d n ii+ 
= ed n u- + (e + l)PiU- + ed n U- + d n u + 

using (|39]l. (|40]) at ft = 0, hence (HOj with 

(42) G=-(e+l)P 1 U- 
Differentiating (|39p twice and using (|36p for u_ and u + we obtain 

o = («_)••- (u+r 

= -a 2 Ag(u_ — u + ) + aQi(d n U- — d rl u + ) + 2a{d n u~ — d n u + ) + ii_ — it+ 

By (|39|) . ([40]) at ft = the first term vanishes and the second equals (e+l)aQi<9 n it_. 
By (|2T))) the third term is (e + l)2a<9„w_ + e2ad n U- — 2aG. Rearranging and using 
(1121) wc obtain J22J) with 



(43) F 2 = -(e + 1)D- 2aed n u- 

D = a(2Pru_ + QxdnU- + 2d n ii. 
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Finally, we differentiate (|40| twice and get, using (|37|) and ([38]) , 

= £'5 n u_ + 2fr(<9 n it_)' + e(<9„u_)" + (<9„u+)" 
= ed n ii- + 2e(Piu_ + Qid n u^ + 9„u_) 

+ P 2 (eii_ + u+) + Q 2 (£<9„u_ + <9„u+) + 2Pi(e?i_ + u + ) + 2<2i(e<9 n w_ + <9„w+) 
+ ed n u- + d n u + 

The terms in the third line of this equation are, using (|39[). (|40p at /i = and (fl"9)) . 

o, 

(e + l)PiU- + + 2Pi((e + l)u_ - F) + 2Q 1 (-ed„u_ + G) 

Using F = — (e + l)a9 n U-, G = — (e + 1)P\U- we obtain 

= ed n u- + 2iA + (e + l)B + ed n u- + d n u+ 
with 

A = Piu_ + d n U- 

B = (P 2 - 2QiPi)u_ + 2P 1 ad n u- + 2P 1 ii^ 
Rearranging, we get (|23[) with 

(44) G 2 = -2eA - (e + 1)P 

Before proceeding we simplify P. Note that 

P 2 -2Q 1 P 1 = -2a 2 HA d + 2a 2 (V g H) ■ V d - 2div g (a 2 W)V d + 4aH (aA d + (V d a) ■ 
= -2dW d (a 2 W)V d + 2div d (a 2 H)V d 
= -2 diva (a 2 Wo)V 9 
where W = W - HI. Thus 

B = -2 diva (a 2 Wo)V a it_ + 2P 1 (ad n u^ + u_) 
We have proved (fT3")) . In order to calculate £, e, we use the following lemma. 

Lemma 12 (Compatibility condition). Lets be aplasmonic eigenvalue with eigenspace 
E. Given functions P, H on dfl, the inhomogeneous system 

Av = ml 3 \ dfl 

V— — v _|_ = P on 9il 

ed n V- + d n v+ = G' on dil 

has a solution v satisfying the decay condition (j4|) iff 

{&,*>-) =e{F,d n w-) 

for all w G E. Here (■, ■) is the scalar product on L 2 (dfl). 

Proof. We reformulate this in terms of Dirichlet-Neumann operators, compare 
Proposition [5] Let g = v+ at dQ. Then the second equation says V- = F + g, 
and the first and third that eAf-(F + g) +N+g = G' . Therefore, a solution v exists 
if and only if the equation 

(eA/1 +N+)g = G' - eN-F 

has a solution g. Since the operator eAf- + A/+ is selfadjoint, this is equivalent to 
(G' — eAf-F, W-) = for all w G E. This is equivalent to the claim. □ 
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We apply the lemma with v = u and w = u and F in f|41 [) . G' = —id n U- + G. 
Using the assumption (u_,9 n u_) = ||u_|| 2 _ = 1 we get 

e = (G,u-) - e(F,d n u-) 

Integrating by parts (G,U-) = (e + l)(divg aVau_, u_) = — (e + l)(aVaw_, Vau_) 
we obtain (|TT]) . 

Now we apply the lemma with v = u and w = u and obtain 

£ = (G 2 ,u-) - e(F 2 ,d n u-) 

Now we sort by like terms. For example, the terms involving u_,it_ are 

-2e(Piw_,u_) - (e + l)(-2)(div a a 2 W V a u-, U-) 

where P\ = — divgaVg. Integrating by parts yields the first line in (|12[) . and the 
other terms are obtained in a similar way. 

3.3. Shape derivative of Dirichlet-Neumann operator. As a side result of 
our calculations above we can calculate the derivative of the Dirichlet-Neumann 
operator Af := Af- when changing the domain as in (J9j> . Since the domain <90 h 
varies, we have to state precisely what we mean by this. Recall from Section 13.11 
the diffeomorphism 2), : — > Q, h , x i-> x + hV(x) where V = an on <90. This 
induces an isomorphism : C°°(dO h ) -> C°°(<90), g H> T^g = g o T h . Then we 
consider 

M h := T*oMn h o (T*)" 1 : C°°(dO) -> C°°(50) 

where A/fi h : C°°(dtlh) — > C°°(dQ,h) is the Dirichlet-Neumann operator for £lh- 
The operator Mh is the same as the one introduced in Section 13.11 Explicitly, it is 
given as follows: Let g S C°°(<90). Let it(-, ft) be the solution of 

(45) Au(-,ft) = in O h 

(46) u{T h (x),h) = g(x) for all x G <90 

(47) and then set (J\f h g)(x) = {d n u){T h {x), h). 

Theorem 13. The derivative at h = of the familiy of Dirichlet-Neumann oper- 
ators defined above is given by 

(48) Af=- div a aV + 2aHAf - AfaAf 

Proof. Differentiate (|35]), (gSJ), (g7J at ft = and use (J37J) (where g h = T h (p)) for 
the right side of (|4"T|) . This gives 

Au = in 

ad n u + u = at 30 

A/g = — diva aV qu + 2aHd n u + d n u 

Rewrite the second equation as u = —ad n u at 90, then together with the first 
equation this gives d n it = Af(—ad n u) = —Af(aAfg), and then the third equation 
gives the claim. □ 

For the exterior Dirichlet-Neumann operator A/+ one has the same formula, 
replacing Af by Af+ everywhere. Same proof. 

Using the second derivative formulas in Section 13.21 one could also derive a for- 
mula for the second derivative of Afh in ft. 
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Note that the principal symbol of the first and third term in (|48| is a|£| 2 , hence 
Af is a pscudodiffcrcntial operator of order at most one, which is to be expected 
since each Afh is a pseudodiffcrential operator of order one. 

3.4. Splitting at first order. We now prove the last statement of Theorem |4j 
which gives also another (though similar) proof of the formula for e. We use the 
characterization in terms of Dirichlet-Neumann operators, Proposition El where we 
use the operators Afh transplanted from tth to f2 as in the previous section. If 
g{x,h) = u(T h (x),h) for x G dil then (e(h)Af-, h + Af +t h)g(h) = for each h. 
Differentiating in h at h = gives 

(sAf- + eAf_ + Af+)g + {sAf- + Af+)g = 

We take the scalar product of this with any g' G Eg :~ {u\g n : u G E}. Using 
sclfadjointncss of Af± and (eAf- + Af+)g' — we get 

((eAf- + eAf- + Af+)g, g') = Vg' G Eg 

or equivalently 

-((eAf-+Af + )g,g')=i(g,g')^ 
A simple calculation using (|48)) for Af = Af± reveals that the left hand side equals 
qi(u, «'), where u' is the harmonic function on fi with boundary values g' , and the 
claim follows. 

4. Examples 

4.1. Planar disk. For the disk f2 = {x G M 2 : \x\ < 1} the only plasmonic 
eigenvalue is e = 1, and any harmonic function in SI whose value and normal 
derivative at dfl are defined (e.g. any function in H s (Sl) with s > 3/2) is a plasmon. 
This follows from the fact that if it- is harmonic in f2 then its Kelvin transform 
it+(x) = u-(j^rs) is harmonic in \x\ > 1. Clearly, u_ and u + have the same 
boundary values at \x\ = 1 and their normal derivatives add to zero, so e = 1. 
The orthogonality in Theorem Q] then shows that there can be no other solutions. 
Choosing u with u\gn non-smooth we see that the regularity statement in Theorem 
[2] does not hold. 

4.2. Perturbed ball. For the unit ball, tt = {x G M 3 : \x\ < 1}, the solutions of 
the plasmonic eigenvalue problem are given by 

k + 1 

u_ fc (r,w) = r k Y k (uj), u +! k(r,Lu) = r^" 1 Y k (u)), e k = 

for k = 0, 1, ... . Here r > 0, w with = 1 are polar coordinates (i.e. x = ru, 
r = \x\, u> = x/\x\). The function Y/. varies over the space of spherical harmonics of 
degree k, i.e. r k Yk(ui) is a fc-homogeneous harmonic polynomial on R 3 . Therefore, 
£fc has multiplicity (2k + 1). 

It is clear that these functions are solutions, and since the Y& span L 2 (S 2 ), there 
can be no further solutions. 

Since the problem is scale invariant, one gets the same e for spheres of arbitrary 
radius. This implies that for a = 1 we should get e k = 0. As a check on our 
calculations we can see this directly: 

From —AgY k = k(k + l)Y k we have for u_ = «_ k 

[ \Vgu_\ 2 = -( ii-Agii- = k(k + 1) / \Y k \ 2 . 
JdQ. Jon Jon 



THE PLASMONIC EIGENVALUE PROBLEM 



21 



Also, d n ii- = d r u- = kYk at r = 1, so the right hand side of (|TT|l equals for a = 1 

[k(k + 1) + e k k 2 ] [ \Y k \ 2 = 0. 
Jan 

For arbitrary a we calculate e in case k = 1 and U-(x,y,z) = z, where £ = 2. 
This choice of u_ (out of the three-dimensional space, spanned by the functions 
x, y, z, of solutions for k = 1) corresponds to an exterior electric field pointing in the 
z-direction, in the physical problem leading to ([J)-©. We parametrize the sphere 
as graph z = ±yl — x 2 — y 2 . Since u_ is 1-homogeneous, we have d r U- = u- = z 
at r = 1. Then |Vait_| 2 = |Vu_| 2 - \d r U-\ 2 = l-z 2 . So the bracket in equals 
(1 — z 2 ) — 2z 2 = 1 — 3z 2 = 3(x 2 + y 2 ) — 2. The volume clement of the sphere is 
dxdyj 'a/1 — x 2 — y 2 . To normalize u we divide by ||it||i = J n |Vu_| 2 = J n 1 = 
vol(fi) = |7r. Finally, e + 1 = 3. Therefore, we get 

9 f „, N 3(x 2 + y 2 )-2 , , 

£i = — / A(x,y)-^=?=J = dxdy, 

^ Jx 2 +y 2 <l yJl-X 2 -y 2 

A(x, y) := o(x, y, V 1 - 2:2 - Z/ 2 ) + ofa, - V 1 - x 2 - y 2 ). 

4.3. Example: Half space. We always assumed ft to be compact. However, a 
similar problem can also be formulated for non-compact Q. For example, consider 
the half space 

il = {{x',x n ) : x' e K n_1 , a-„ < 0} 

We then consider the problem ([I|)-((3]), but we need to replace the decay condition 
(j4|. In view of the regularity discussion it is natural to consider H 1 / 2 boundary 
data, hence we demand a uniform H 1 / 2 condition: there is a constant C so that 

(49) / \u(x',x n )\ 2 dx' + / \u(x',x n )\ \V x 'u(x', x n )\dx' < C for all x n 

JR"-1 it"" 1 

In particular, we can take the Fourier transform in the x'-variables, u(£',x n ) = 
/ Rn -i e~ lx u(x',x n )dx', and Ait = yields (d 2 n — |£'| 2 )u = for x n ^ 0, with 
solutions u± = o±e-l £ 'l x " + b±e^'\ x " (if £' ^ 0). Then 0SJ implies 6+ = a_ = 0, 
and ([2]) gives 6_ = a+, so 

u(x',x n ) = J e ix '<'h{0^ WlM ^' 

for a function h on R™ -1 satisfying / Rn _ 1 \h(C)\ 2 (l + \?\)d£' < oo. We have h = 
{2-K)- n g for the boundary value g (x') =u(x',0). Also d n u±(x', 0) = t / e lI ' f '|f 
(distributional Fourier transform), which shows that e = 1 for any g. 

Therefore, e = 1 is the only plasmonic eigenvalue for the half space. Incidentally, 
the explicit formulas show that the Dirichlet-Neumann operator is precisely the 
Fourier multiplier by |£'|, i.e. Af~g = {\£'\gj- 
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